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! ABSTRACT 

ly-j ■ Double neutron star (NS-NS ) binaries, analogous to the well known Hulse-Taylor pulsar 

PSR 1913+16 (Hulse & Taylor 1975), are guaranteed-to-exist sources of high frequency gravi- 
tational radiation detectable by LIGO. There is considerable uncertainty in the estimated rate 
of coalescence of such systems (Phinney 1991; Narayan et al 1991; Kalogera et al 2001), with 
conservative estimates of ~ 1 per million years per galaxy, and optimistic theoretical estimates 
one or more magnitude larger. Formation rates of low-mass black hole-neutron star binaries may 
be higher than those of NS-NS binaries, and may dominate the detectable LIGO signal rate. 
Rate estimates for such binaries are plagued by severe model uncertainties. Recent estimates 
suggest that BH-BH binaries do not coalesce at significant rates despite being formed at high 
rates (Portegies Zwart & Yungelson 1998; De Donder & Vanbeveren 1998). 

We estimate the enhanced coalescence rate for BH-BH binaries due to weak asymmetric kicks 
during the formation of low mass black holes like Nova Sco (Brandt, Podsiadlowski & Sigurdsson 
1995), and find they may contribute significantly to the LIGO signal rate, possibly dominating 
. £h ! the phase I detectable signals if the range of BH masses for which there is significant kick is 

broad enough. For a standard Salpeter IMF, assuming mild natal kicks, we project that the R§ 
merger rate (the rate of mergers per million years in a Milky Way-like galaxy) of BH-BH systems 
is ~ 0.5, smaller than that of NS-NS systems. However, the higher chirp mass of these systems 
produces a signal nearly four times greater, on average, with a commensurate increase in search 
volume. Hence, our claim that BH-BH mergers (and, to a lesser extent, BH-NS coalescence) 
should comprise a significant fraction of the signal seen by LIGO. 

The BH-BH coalescence channel considered here also predicts that a substantial fraction of 
BH-BH systems should have at least one component with near- maximal spin (a/M ~ 1). This 
is from the spin-up provided by the fallback material after a supernova. If no mass transfer 
occurs between the two supernovae, both components could be spinning rapidly. The waveforms 
produced by the coalescence of such a system should produce a clear spin signature, so this 
hypothesis could be directly tested by LIGO. 

Subject headings: black hole physics — binaries: close — stars — gravitational radiation. 
1. Introduction 
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field of gravitational radiation astronomy is set to 
prosper, if there are detectable sources in the lo- 
cal universe. A considerable amount of effort has 
been directed towards identifying potential grav- 
itational radiation sources, and the relative rate 
of contribution of these sources to the anticipated 
signals. The canonical scenario envisioned is the 
final coalescence of a binary neutron star system 
(NS-NS), where the pair's orbital energy has been 
radiated away by gravitational waves. 

NS-NS merger can produce a copious signal for 
LIGO, if they occur at high enough a rate locally 
- current estimates for LIGO phase I place the 
maximum detection radius for such an event at 
~20Mpc or less (Kalogera et al 2001). The merger 
rate of NS-NS binaries can be estimated from 
observed systems (Phinney 1991; Narayan et al 
1991), such estimates are plagued by small number 
statistics and possible observer biases (Kalogera 
1998; Kalogera et al 2001). Alternatively, the rate 
can be estimated from ab initio theoretical models 
(see review by Grishchuk et al 2001, also Portegies 
Zwart & Yungelson 1998; Belczynski et al 2001; 
Fryer et al 1999; Bloom, Sigurdsson & Pols 1999). 

A number of authors have explored binary pop- 
ulation synthesis models, making a number of 
assumptions about the input physics, leading to 
rates consistent with observational constraints, 
but uncertain by 1-2 orders of magnitude (Porte- 
gies Zwart & Yungelson 1998; Kalogera 1998; 
Kalogera & Lorimcr 2000; Bloom, Sigurdsson 
& Pols 1999; De Dondcr & Vanbeveren 1998; 
Belczynski & Bulik 1999; Lipunov, Postnov & 
Prokhorov 1997; Brandt & Podsiadlowski 1995; 
Belczynski et al 2001; Fryer et al 1999). Con- 
ventionally, the rates are expressed in terms of 
i?6, the merger rate per million years per Milky 
Way-like galaxy, assuming normal rates of star 
formation. Then the rate in the local universe 
is the integrated rate over the number density of 
galaxies, the rate per galaxy, scaled to the Milky 
Way rate, per detection volume. Kalogera et al 
(2001) find a LIGO I event rate of 3 x 10~ 4 y" 1 
for NS-NS mergers within 20 Mpc, assuming 
Rq = 1. For a given detector sensitivity, higher 
mass coalescences are detectable to larger vol- 
umes, with the detection distance, oc M^f , 

where M chlrp = (mim 2 ) 3 '' 5 /(mi +m 2 ) 1/b (Thorne 
1994). Since black holes are expected to have 
masses several times larger than neutron stars, 



and the event rate scales as d\ (for dh small com- 
pared to the size of the universe), a black hole 
coalescence rate of order R 6 implies event rates 
2-3 orders of magnitude higher. 

The problem of accurately estimating the coa- 
lescence rate of compact binaries can be appreci- 
ated by noting that the type II supernovae rate 
in the Milky Way is 1-2 per century, implying 
that the rate of type II supernova in binaries is 
about 10~ 2 y _1 , given a 50% binarity rate. Esti- 
mates of the coalescence rate are canonically close 
to i? 6 = 10~ 6 y" 1 per Milky Way. So the branch- 
ing ratio for type II supernovae to form merging 
systems is of the order 10~ 4 . Calculating the mean 
rate, averaged over all scenarios for coalescence, 
in the local universe, is hard, with small errors of 
assumption about the physics of stellar evolution 
leading to large fractional changes in the branch- 
ing ratios of particular channels for mergers. 

Theoretical models require a series of assump- 
tions, about the mass function of high mass stars, 
the cut-off points for which zero age masses lead 
to NS or BH formation, the binary fraction and 
mass ratio distribution, the amount of mass loss 
during (binary) stellar evolution, and the ampli- 
tude and distribution of natal kicks. Secondary 
assumptions, that are usually not explored in de- 
tail, include the dependence of all of the above on 
metallicity and environment in which the massive 
stars form; and the possibility that the NS/BH 
cut-off, and the natal kick, may depend on stellar 
rotation and magnetic fields. In fact, it is not im- 
plausible that the NS/BH formation boundary is 
"fuzzy", that there is not a sharp border in zero- 
age mass between stars that form neutron stars 
and those which form black holes. A further po- 
tential confounding effect, is that natal kicks are 
generally assumed to be random, but may in prin- 
ciple be correlated with some macroscopic prop- 
erty of the progenitor star, such as rotation (e.g. 
Spruit & Phinney 1998; Pfahl et al 2001). 

There arc two arguments against a high rate 
of coalescence of compact binaries: most high 
mass binaries become unbound during the near 
instantaneous mass-loss at supernova; the second 
problem is that of the 'kicks' believed to be in- 
flicted upon the NS from an asymmetrical super- 
nova blast. While the range of energies imparted 
by the kick is a subject of intense debate (Lyne 
& Lorimer 1994; Hansen & Phinney 1997; Tau- 
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ris et al 1999; Cordes & Chernoff 1998; Fryer et al 
1998; Arzoumanian, Chernoff & Cordes 2001), the 
resulting velocity change is believed to be of the 
order 250 km sec -1 for a 1.4 M Q neutron star. A 
properly aligned kick can allow the system to re- 
main bound, albeit with dramatically altered or- 
bital parameters, but a random kick is unlikely to 
be directed so fortuitously, with the extra energy 
making a bound final state even less likely. 

An alternative to the NS-NS scenario is one in 
which both stars are of sufficient mass to end up 
as black holes. The minimum mass required for 
black hole formation is not well known, and is 
heavily influenced by the evolutionary history of 
the black hole candidate (see Fryer 1999; Fryer & 
Kalogera 2001). Accretion induced spin-up of the 
star prior to collapse could result in a higher mini- 
mum mass threshold, while the presence of strong 
winds could drop the progenitor mass below the 
required minimum. Finally, unless an unusually 
flat IMF is assumed (as apparently seen in cer- 
tain starburst galaxies, see, for example, Doane & 
Mathews 1993), few of these objects will be pro- 
duced in the first place. 

Nevertheless, precursors to black hole binary 
systems are more likely to survive both super- 
novae due to the smaller fractional mass loss dur- 
ing the event. Depending upon the mass of the 
star, a significant amount of material can fall 
back upon it just after the blast. Indeed, a suf- 
ficiently massive progenitor may not undergo a 
supernova at all, reaching a point where nearly 
all of the material is quickly reabsorbed (Fryer 
1999). After losing a significant fraction of mass 
to stellar winds and (possibly) two supernovae, a 
typical bound binary system will have expanded 
considerably, with most P or b > 10 days. Long- 
period systems do not decay in less than a Hub- 
ble time from gravitational radiation, and hence 
will never reach a frequency range useful to LIGO. 
However, including asymmetric supernova kicks 
changes the picture significantly. Assuming that 
the total momentum change is approximately the 
same for all supernovae, the velocity imparted to 
the nascent black hole will be scaled to its mass, 
via Avbh — AvnsMns/Mbh (see Grishchuk et 
al 2001 for review). For a typical 7M© hole, this 
gives a kick of the order ~ 50 kms -1 , compara- 
ble to or larger than the orbital speed, but not so 
large as to always rip the system apart. Systems 



that remain bound are generally in highly eccen- 
tric orbits, which expedites their merger through 
emission of gravitational radiation. 

Ultimately, the point is this: the trace of 
the metric perturbation from gravitational waves, 
varies as h cx M^[f rp . As the maximum detection 
radius scales linearly with h, the maximum search 
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volume scales as V cx M c ^ irp . The implied sev- 
eral hundred-fold increase in search volume means 
that, even if only a small fraction of these systems 
can merge in t < 1/Hq, they may dominate the 
detected signal. 

Our approach is to concentrate on estimating 
the event rate from one particular coalescence 
channel. The total branching ratio for coales- 
cence is of course the sum of all possible channels, 
and it is possible that other channels contribute 
significantly to the event rate, possibly dominat- 
ing the total rate; for example, dynamical evolu- 
tion of cluster binaries may create new channels 
for merger with high coalescence rates (Sigurds- 
son & Hernquist 1993; Portegies Zwart & McMil- 
lan 2000), or accretion induced collapse of neutron 
stars with soft cores may lead to enhanced rates 
(Bethe & Brown 1998). Here, rather than trying 
to estimate the total event rate, we make an obser- 
vationally motivated estimate of the rate for one 
particular channel. 

2. Example of Nova Sco 

Nova Sco 1994 (GRO J1655-40) is a strong can- 
didate for being a black hole binary (Bailyn et al 
1995; Orosz et al 1997), with the primary being a 

6.3 ± 0.5 M black hole (Greene, Bailyn & Orosz 
2001; Shahbaz et al 1999). Perhaps its most re- 
markable feature is an unusually high space ve- 
locity, whose lower limit of 106 km s _1 was, un- 
til recently, several times greater than any other 
known black hole transient (Shahbaz et al 1999; 
Brandt, Podsiadlowski & Sigurdsson 1995). The 
likely true space velocity is greater by a factor of 
•s/3, adjusting for mean projection effects. 

A number of scenarios have been put forward 
to explain the unusual speed of Nova Sco. Brandt, 
Podsiadlowski & Sigurdsson (1995) point out that 
the momentum component of Nova Sco along the 
line of sight is comparable to that of a single 

1.4 M neutron star, having received a natal kick 
in the range 300-700 km s _1 . This is not an un- 
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reasonable value for a neutron star kick (Lyne & 
Lorimer 1994; Cordes & Chernoff 1998; Fryer et al 
1998), and so lends strength to the possibility that 
Nova Sco can be explained by the primary experi- 
encing a natal kick prior to formation of the black 
hole. Invoking Blaauw-Boersma kicks (Blaauw 
1961) as the sole acceleration mechanism invites 
difficulty as the low mass of the secondary means 
that most Blaauw-Boersma kicks strong enough 
to give the observed speed would also disrupt the 
binary system completely (Brandt, Podsiadlowski 
& Sigurdsson 1995). This is not to say that such 
a scenario is impossible (see, for example, Nele- 
mans ct al 1999). However, the space velocity 
measurement of the X-ray nova XTE J1118+480 
has provided a further example of a black hole 
binary with a high- velocity Galactic-halo orbit. 
Mirabel et al. (2001) used the VLBA to obtain 
a precise proper motion for the system, calculat- 
ing a speed of 145±35 km s -1 , with respect to 
the local standard of rest. The system primary 
has a mass function of 6.0±0.4M Q (McClintock et 
al 2001), with a faint (~ 19 mag) optical counter- 
part. To accelerate this system to the observed pe- 
culiar velocity using Blaauw-Boersma kicks alone, 
roughly 40 M Q of material would have to be ex- 
pelled during stellar collapse, an implausibly large 
amount of matter. This issue remains contentious, 
however, as calculations in Nelemans et al (1999) 
show that the velocity of Nova Sco, at least, might 
be explained by a symmetric supernova where the 
black hole progenitor lost more than half of its 
mass. This is marginally possible, but becomes 
much more tenuous if Nova Sco turns out to have 
a significant transverse velocity. 

One principal difficulty with the natal kick 
premise is arranging for sufficient neutrinos to 
drive a supernova explosion prior to being trapped 
by the formation of an event horizon, ejecting neu- 
trinos or other material asymmetrically. Few neu- 
trinos escape if the horizon forms over the dynam- 
ical time of the collapsing core (Gourgoulhon & 
Haensel 1993). The drop in the resulting neutrino 
heating of the envelope allows most or all of the 
envelope to fall back onto the nascent black hole 
(see, for example, the hydrodynamical simulations 
of Janka & Mueller 1996) , hence preventing the su- 
pernova. However, there is now strong evidence 
that a supernova must have taken place in the 
Nova Sco system, from estimates of metal abun- 



dances in the atmosphere of the secondary star. 
Israelian et al (1999) show that the secondary, an 
F3-F8 IV/III star, possesses a dramatic enhance- 
ment (factors of six to ten) in a elements such 
as oxygen, magnesium, silicon and sulfur (and see 
also the discussion in Podsiadlowski et al 2001). 
Interestingly, no significant enhancement in iron 
was found. The implication is that the progenitor 
of the black hole did experience a supernova dur- 
ing its formation, and a significant amount of the 
former atmosphere was captured by its compan- 
ion. As most of the iron core went on to collapse 
to a singularity, no enhancement of this particular 
element is seen. The question is then how long 
after the supernova it took for the black hole to 
form. One possibility is that the star, at the in- 
stant of the supernova, experienced sufficient ro- 
tational support to avoid collapse for an extended 
period of time, until spin-down allowed the system 
to collapse to a black hole some time later. An- 
other possible scenario has some of the material 
cast off in the supernova explosion recaptured by 
the newly-formed neutron star, elevating its mass 
above the threshold for collapse. In cither case?, a 
supernova is seen, with a corresponding asymmet- 
ric natal kick. More massive progenitors would 
form a horizon directly, with no possibility for a 
kick. We discuss limits for this behaviour in the 
next section. Nova Sco 1994 and XTE J1118+480 
provide strong evidence that black holes do, in 
fact, experience natal asymmetric kicks, at least 
under some conditions. The resulting orbital ec- 
centricities tend to accelerate binary coalescence 
from the enhanced radiation emitted at every pe- 
riastron passage. 

3. Population synthesis 

We make use of a binary evolution code devel- 
oped by (Pols & Marinus 1994), and modified for 
use in NS-NS systems by (Bloom, Sigurdsson & 
Pols 1999). Our extension of the code allows for 
evolution to the black hole state, with assump- 
tions about the mass function of such objects at 
the time of collapse. 

Initially, the code chooses the mass of the pri- 
mary from a given mass function. For this work, 
we have chosen two power law IMFs, with indices 
a of —2.0 and —2.35 (the latter, of course, be- 
ing the Salpeter IMF). In both cases, we estab- 
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lished a lower cutoff of 4M Q , confining the code to 
an interesting range of initial masses; i. e., where 
at least one supernova is possible. Our primary 
stellar models are taken from Maeder & Meynet 
(1989). The helium star models used are a mix of 
models from Habets (1986) and Paczynski (1971). 
We assume that the mass ratio distribution be- 
tween the two components is flat. After choos- 
ing an initial separation a and eccentricity e after 
Pols & Marinus (1994), the code evolves each bi- 
nary system until both components have reached 
their final degenerate form, accounting for mass- 
transfer-induced stellar regeneration and stellar 
winds. During the common-envelope phase, the 
orbit is circularised, and the orbital energy is re- 
duced by the binding energy of the envelope di- 
vided by the common-envelope efficiency parame- 
ter, which we take to be 0.5. In other words, the 
orbital energy is reduced by twice the envelope 
binding energy Neutron stars are formed from 
progenitors with ZAMS masses of between 8 and 
20 Mq, inclusive, and are always given a mass 
of 1.4 Mq. More massive stars end up as black 
holes. As noted in the Introduction, this bound- 
ary is likely to be "fuzzy", i.e., not a monotonic 
function of the progenitor mass, as it is strongly 
coupled to the spin state of the star prior to col- 
lapse, a quantity which our code simply does not 
track. Even assuming this was known to perfect 
accuracy, it is far from trivial to estimate the ef- 
fects of magnetic field and rotational support vis- 
a-vis the compact object's end state. 

The black hole mass function (i. e., the post- 
collapse mass of a BH, given its mass just prior to 
the explosion) is highly speculative at this point, 
and is almost certainly not merely a function of 
initial mass, but also of angular momentum, to 
the extent that this determines the fraction of ma- 
terial falling back onto the collapsing star. In or- 
der to experience a kick, the black holes formation 
must be delayed somewhat, either from rotational 
support, or because event horizon formation oc- 
curs only after delayed fallback of mass initially 
ejected from the core. Fryer (1999) has performed 
core-collapse simulations in order to explore the 
critical mass for black hole formation, and their 
final masses. As a best working guess, we have 
constructed a mass relation from a quadratic fit 
to the limited data set found in Fryer (1999). The 
mass of the black hole at formation (Mbh) is re- 



lated to the ZAMS mass of the progenitor (Mo) by 
M B h = (Mo/25) 2 x 5.2M Q . We have also chosen 
a simple criterion for whether a black hole will re- 
ceive an asymmetric neutrino kick during collapse; 
namely, all objects below 40 M Q (referring to the 
ZAMS mass) experience a random kick. Above 
this limit, objects collapse directly to a black hole, 
with no kick. We explore the effect of ignoring 
black hole kicks in the next section. 

As a last step, the results of the code are nor- 
malised to the supernova rate of the Milky Way, 
0.01 yr _1 . Here, the code produces some arbi- 
trarily large number of complete evolutionary se- 
quences, including a large number of supernova 
events. Scaling all these events to the Milky Way 
SNR allows us to generate event rates out of data 
not previously ordered in time. This effectively 
simulates a constant star-formation rate, where 
the population of merger candidates has reached 
a dynamic equilibrium. 

4. Results 

Tables 1 and 2 summarize the investigation. 
For each IMF, six simulations (generating 10 6 bi- 
nary systems each) were run, varying two param- 
eters: if max, the upper limit of the imparted kick 
speed (acting on a 1.4 M© NS), and <r, the dis- 
persion in the Maxwellian distribution of the kick 
speed. Maximum kick values were chosen at 500 
km s _1 (for a = 90 km s _1 and a — 190 km s _1 ), 
and 1000 km s" 1 (for a = 450 km s _1 ). The inclu- 
sion of the 90 km s _1 kick intensity is based upon 
recent observational work on the velocity distribu- 
tion of pulsars by Pfahl et al (2001), which sug- 
gests a bi-modal kick distribution centered around 
<~ 100 and ~ 500 km s _1 . Further discussions 
concerning this bi-modality can be found in Ar- 
zoumanian, Chernoff & Cordes (2001); Cordes & 
Chernoff (1998); Fryer et al (1998). This range has 
been touched on in previous work (Grishchuk et 
al 2001), but we have looked at it specifically as a 
new observationally-motivated choice in kick mag- 
nitude. The following columns show the Rq rate 
of BH-BH , BH-NS , NS-NS and BH-PSR pairs 
generated in the sample. Rq is a rate of 1 event 
per Myr per Milky Way galaxy, found by scaling 
event rates to the estimated local supernova rate 
of 0.01 yr -1 . One must take into account that not 
all stars (and hence not all supernovae) occur in bi- 
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naries; hence a binarity fraction must be assumed, 
and we take that value to be 0.5 throughout (i. e. 
roughly two-thirds of supernovae occur in a binary 
system). The alternating columns show the R§ of 
such pairs expected to coalesce through emission 
of gravitational radiation in less than the Hubble 
time (here, taken as 10 Gyr). Table 1 details the 
formation and merger rates when black holes and 
neutron stars are given kicks as described above. 
Table 2 gives the relevant rates when black holes 
are not allowed to have asymmetric kicks of any 
kind (even if below 40 M Q ). 

Figure 1 show examples of the generated data 
sets. Here, the separations of the BH-BH and BH- 
NS systems after the second supernova (i. e., stel- 
lar evolution has ended) are plotted against the 
post-mortem orbital eccentricity. Dots represent 
systems stable against orbital decay, whilst trian- 
gles represent those systems which will undergo 
coalescence in less than a Hubble time. Figure 2 
shows the distribution of masses in surviving BH- 
BH binaries, irrespective of whether the system is 
unstable to orbital decay.. A clear bifurcation in 
mass is seen, with the vast majority at around 
5-6 Mq, and another, smaller concentration at 
roughly 35 M Q . The low end is comprised of black 
holes with 20-40 M ZAMS progenitors, which 
lose much of their mass through a supernova (and 
hence experience a natal kick). The higher mass 
group are from more massive progenitors, and col- 
lapsed to a black hole without a supernova or the 
resulting kick. The histogram in figure 3 shows 
the "chirp" mass distribution for all BH-BH sys- 
tems that eventually merge from gravitational ra- 
diation. As the wave amplitude scales like M^lf , 
this plot, coupled with the sensitivity of the LIGO 
detector, allows us to estimate a detection rate for 
BH-BH binaries. By way of comparison, two 1.4 
Mq objects produce a chirp mass of 1.22 M©. 

As can be seen from column 8 in tables 1 and 2, 
most NS-NS pairs that remain intact coalesce in 
r < 1/Ho, irrespective of IMF and kick strength. 
This is perhaps not surprising, since NS-NS sys- 
tems must generally be tightly bound to remain 
intact through the second supernova. For NS- 
NS systems, we estimate the Galactic merger rate 
to range from ^0.04-8 Myr -1 , with a strong de- 
pendence on the severity of the natal kick. Given 
a likely LIGO I detection radius of <~ 20 Mpc for 
such systems, this would imply a conservative de- 



tection rate of 0.01-3 xlO^yr" 1 (Phinney 1991; 
Kalogera & Lorimer 2000), with a likely mean rate 
of ^10~ 3 yr -1 . These estimates are an order of 
magnitude lower than those stated in Grishchuk 
ct al (2001). 

A more unexpected result was the relative fre- 
quency of BH-NS pairs relative to bound NS sys- 
tems. While far fewer BH-NS systems are formed, 
their mass allows them to weather two supernovae 
more frequently than two neutron stars. From an 
estimate of the BH-NS formation rate, we can say 
something about the expected frequency of BH- 
PSR systems. For purposes of discussion, we di- 
vide these systems into three categories, based on 
formation mechanism. The first type results from 
the standard scenario: a massive black hole forms 
first, with a regular short-lived pulsar at a later 
epoch. The second possibility is that the black 
hole progenitor transfers a substantial amount of 
matter on to the neutron star progenitor early in 
the evolution of the system. This results in the 
neutron star forming first, and allows the black 
hole progenitor to continue to spill matter through 
its Roche lobe on to the slowing pulsar, thereby 
recycling it. The last possibility is that the sys- 
tem from scenario two is disrupted by the black 
hole formation, resulting in a (possibly recycled) 
pulsar and single black hole. This channel is not 
relevant for our purposes, as it ceases to be a po- 
tential LIGO source. The Rq rates for the remain- 
ing formation channels are shown in tables 1 and 
2. The BH-NS column covers all formation chan- 
nels for a BH-NS system, whereas the BH-PSR2 
column considers only that subset of BH-NS sys- 
tems where a recycled pulsar is formed first, via 
channel two, as described above. Surprisingly, of 
the systems that go on to form a BH-NS pair, 
around 40% experience the reversal described in 
case two above, where the initially more massive 
star forms a compact object last. Only 0.5% of 
the BH-BH progenitor systems experience a simi- 
lar reversal. In the first case, an ordinary, short- 
lived (~ 10 6 yr) pulsar is created. In the latter, 
a fast, long-lived (~ 10 9 yr) millisecond pulsar 
should often be the result. It is likely that type 
II BH-PSR systems with recycled pulsars would 
be observable as radio objects. The small pre- 
dicted formation rate of such objects, however, is 
consistent with having not yet been seen among 
known pulsars. However, their formation rate, 
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while 30-60 times smaller than that of "type 1" 
BH-NS systems, should mean that they are con- 
siderably more common than black holes with nor- 
mal pulsar companions, being longer-lived. Also, 
we note that the formation rate of "type 2" BH- 
PSR systems is more severely affected by kick 
intensity, as strong kicks may widen the system 
enough to prevent mass transfer on to the neu- 
tron star, leaving an ordinary BH-PSR system. 
We conclude that the galactic scale height for such 
systems should be substantially smaller than that 
of ordinary BH-NS and BH-PSR systems. Further 
details of the BH-PSR systems (including a predic- 
tion of orbital distributions) will be deferred to a 
future work (Sipior & Sigurdsson, in prep). Note 
finally that the chirp mass of a typical coalesc- 
ing BH-PSR system is roughly twice as large as a 
NS-NS binary. As d L scales with M c 5 ^ 6 rp , this im- 
plies at least a five-fold increase in search volume. 
Adjusting the estimated Rq rates for this implies 
that BH-NS coalescence should be seen with sub- 
stantially greater frequency than NS-NS binaries 
at all kick intensities, with a conservative maxi- 
mum of 7 x 10~ 3 LIGO I detections per year for 
low kick strengths. A more likely "best bet" rate 
is on the order of 2-3xl0~ 3 yr -1 , with a mix of 
high (a = 450 km s _1 ) and low (a = 90 km s _1 ) 
kick velocities, with 60% of kicks drawn from the 
former, the balance from the latter (as suggested 
by Arzoumanian, Chcrnoff & Cordes 2001; Pfahl 
et al 2001). 

BH-BH systems are shown here to be far more 
common than either NS-NS or BH-NS hybrids. 
Only a small fraction of these will merge in a Hub- 
ble time, as the larger mass permits the system to 
remain bound at higher separations. As well, the 
greater mass of BH-BH systems implies a more 
gradual response to an increasing kick parame- 
ter. Indeed, the Rq rate for bound BH-BH sys- 
tem formation remains high (more than thirty per 
Myr in the Milky Way), even when kicks are per- 
mitted to reach 10 3 km s _1 on a 1.4 M Q object. 
Figure 3 shows the distribution of chirp masses 
for all BH-BH systems that coalesce in less than 
10 Gyr. As can be seen, there is a concentra- 
tion of chirp masses around 4 M , implying a 
twentyfive-fold increase in search volume. Allow- 
ing for this, we anticipate a LIGO I detection 
rate of 0.3-12xl0~ 3 yr" 1 , with a likely rate of 3- 
4x 10~ 3 yr^ 1 , given the two-component kick speed 



distribution described above. Figure 4 shows the 
strong correlation between the final system veloc- 
ity (a function of the kick magnitude and compo- 
nent masses) , and the time required for the system 
to coalesce due to emitted gravitational wave en- 
ergy. A clear dynamic criterion for coalescence 
emerges, with slower systems rarely merging in a 
Hubble time. This is essentially a selection effect, 
as the kick speed must be of a magnitude compa- 
rable to the orbital speed of the system to have a 
dramatic effect on the merger rate. Relative to the 
orbital speed, a weak kick will have no appreciable 
effect, while a strong kick will disrupt the system. 
So, the kick speeds that bring about more rapid 
mergers are necessarily tuned to the orbital speed 
of the system at the second supernova. 

5. Discussion 

We claim a "best guess" LIGO I detection rate 
of 3-4 xl0~ 3 yr _1 for BH-BH coalescence events, 
2-3xl0" 3 yr- 1 for BH-NS events, and ~10~ 3 NS- 
NS coalescence per year. LIGO II, with an an- 
ticipated order of magnitude increase in sensitiv- 
ity, should encompass roughly 10 3 times the search 
volume, with a proportional increase in event de- 
tection rates. Unless LIGO I is quite lucky, it 
seems unlikely that mergers from this channel will 
be detected until the advent of LIGO II in a few 
years. With LIGO II, we can anticipate an event 
rate of many per year. These are essentially consis- 
tent with the most pessimistic detection rates cal- 
culated by Grishchuk et al (2001). There are sig- 
nificant differences between our results and those 
in Portegies Zwart & Yungelson (1998), primarily 
regarding BH-BH coalescence, where a negligible 
merger rate was found. This discrepancy is likely 
due to our lower cutoff for black hole formation, 
where we adopt a value of 20 M Q , against the 40 
M Q cutoff adopted by the latter. While recent 
models of Fryer (1999) have established a lower 
black hole mass cutoff, more work is needed to re- 
solve this area of contention. LIGO-measured co- 
alescence rates should provide a clearer resolution 
on this point. 

It is instructive to compare our results with 
those obtained in two recent population synthe- 
sis studies. Fryer et al (1999) performed detailed 
synthesis calculations in the context of black hole- 
accretion driven gamma ray bursts. As a result, 
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BH-BH binaries are not considered; however, BH- 
NS and NS-NS systems are. For the Salpeter IMF, 
and a Gaussian kick distribution peaked at 100 km 
s _1 , FWHM of 50 km s _1 (roughly correspond- 
ing to our o = 90 km s _1 kicks), the Rq rate of 
BH-NS formation is identical (Rq — 16). We pre- 
dict many more NS-NS (Rq of 11, as opposed to 
4.2). This discrepancy arises from the inclusion 
of hypercritical accretion in Fryer et al scenario I 
for NS-NS formation. The result is that all neu- 
tron stars that pass through a common envelope 
phase become black holes, removing a significant 
channel of NS-NS pair formation. At higher kick 
intensities, the BH-NS formation rates remain in 
good agreement, whilst the NS-NS rates become 
more disparate. This is likely because the two high 
speed kick distributions are not the same, and the 
discrepancy affects BH-NS systems less dramati- 
cally because of the larger mass involved. 

Very recently, the paper of Bclczynski et al 
(2001) performed a similar comprehensive study, 
considering a wide variety of kick paradigms and 
population synthesis models. Our anticipated 
LIGO I detection rates are consistently an order 
of magnitude lower than those reported by the au- 
thors for their "standard model" . Though the kick 
distribution chosen for this standard model is more 
intense than in our own work (two Maxwellian dis- 
tributions, with a = 175 km s _1 and a = 700 km 
s _1 , with 80% of kicks drawn from the former, 
and the balance from the latter distribution), the 
mean neutron star mass is allowed to be much 
higher (up to 3 M©, allowing more systems to 
remain bound, and shed gravitational radiation 
more quickly. Thus, the overall effect is to enhance 
the anticipated NS-NS and BH-NS merger rates. 
It should be noted that our calculated merger 
rates fall within the range of merger rates calcu- 
lated by Belczynski et al when all posited model 
assumptions are considered. The difference be- 
tween our results and the standard model of the 
authors for BH-BH systems is larger, by roughly 
two orders of magnitude (though our results still 
fall within the wider range resulting from consid- 
ering all the models presented). Our model as- 
sumptions are greatly discrepant at one point in 
particular; namely, the assumed mass function of 
newly-formed black holes as a function of ZAMS 
mass. Where wc take a rough interpolation be- 
tween model results presented in Fryer (1999), the 



authors use a much more sophisticated relation- 
ship between the ZAMS mass and the initial mass 
of the resulting black hole. This mass function 
tends to produce substantially more massive holes 
than those generated from our interpolation, and 
we believe this to be the principal cause of the 
discrepancy in merger rates for BH-BH binaries. 

After the neutron star/black hole mass cutoff, 
the second pivotal assumption concerns the distri- 
bution of natal kick speeds, and indeed, whether 
black holes can receive kicks under any circum- 
stances. The cases of Nova Sco 1994 and XTE 
J1118+480 show that our natal kick assumptions 
are plausible. However, we suffer from a paucity 
of data concerning the frequency and scaling (with 
respect to mass) of black hole kicks. There is 
mounting evidence that neutron star kick distri- 
butions are bi-modal (Arzoumanian, Chernoff & 
Cordes 2001), allowing many more NS-NS systems 
to survive and merge after the second supernova. 
The greater mass of black holes means that uncer- 
tainties in these distributions have less effect on 
the coalescence of BH-NS and BH-BH systems, as 
can be seen in table 1. Coupling this fact with 
the "residual" merger rate one gets when black 
holes are not allowed natal kicks at all (cf. table 
2), we argue that the effects of uncertainty in the 
natal kick models will not affect the predicted Rq 
rates of BH-BH and BH-NS mergers by more than 
factors of a few. Interestingly, the large number of 
BH-BH binaries which remain bound suggests that 
another merger channel; namely, a merger due to 
hardening resulting from external dynamical per- 
turbations (Sigurdsson & Hernquist 1993). We do 
not quantify this effect here, but note that the 
extent and duration of a star-formation episode 
could have a notable effect on the rate at which 
such interactions occur. 

Black holes experience natal kicks only when 
formation of the event horizon is delayed, from 
fallback material and/or initial spin support. 
Therefore it is reasonable to posit a correlation 
between natal kicks and a high initial spin for the 
nascent black hole. In a BH-BH system, the first 
hole is gradually spun down by interaction with 
the stellar wind of its companion, and will likely 
be a slow rotator even if it began with maximal 
spin. The second hole will not experience this, 
and should retain its high spin. We therefore pre- 
dict that, for BH systems that experience natal 
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kicks, LIGO will see a slow rotating primary, with 
a (generally less massive) black hole companion 
with near maximal spin. If the components do 
not experience natal kicks, we would anticipate 
low spins all around, with the coalescence rate 
primarily determined by the wind mass loss rate 
and the boundary criterion for black hole forma- 
tion. This spin signature is a directly testable 
hypothesis — LIGO II should provide a large array 
of high-quality sources to verify this prediction, 
and allow for careful testing of gravitational wave 
phenomenology. 

A similar spin correlation may be found in BH- 
PSR systems, with the two categories described 
above having different anticipated signatures. The 
type I systems (black hole first, followed by stan- 
dard pulsar) should have two slow rotators, as no 
pulsar recycling is possible. Type II systems (neu- 
tron star forms first, is spun up by black hole pro- 
genitor, black hole forms) should have at least one 
rapid rotator, and possibly two, if the black hole 
is spun up during its formation, as above. These 
spin correlations may be demonstrated with high- 
S/N data from LIGO II. More will be said about 
the spin signatures of BH-PSR systems in an up- 
coming paper (Sipior & Sigurdsson, in prep), in 
addition to a study of the orbital parameters of 
such systems, and the ratio of BH-PSR to NS-PSR 
systems, which is relevant for new tests of General 
Relativity, should these systems be detected. 

Model uncertainties, especially involving the 
details of black hole formation, still dominate 
the estimates of coalescence rates to be seen by 
LIGO. Nonetheless, it is becoming clear that BH- 
BH mergers form a significant channel, compara- 
ble to the rates anticipated for NS-NS systems. 
While LIGO I detection rates appear too small 
to catch any merger events from this channel, we 
may still be lucky enough (at the level of a few 
percent) to serendipitously catch events from this 
channel (especially given incremental upgrades in 
the sensitivity of LIGO I throughout its opera- 
tional lifetime). LIGO II should easily see a large 
number of these events every year. 
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Table 1 



a 


a (km s 1 ) 


BH-BH 


T <T h 


BH-NS 


T <T h 


NS-NS 


T <T h 


BH-PSR 2 


T < T h 


-2.0 


90 


80 


1.6 


26 


8.2 


11 


8.4 


0.62 


0.44 


-2.0 


190 


68 


0.59 


5.9 


2.5 


1.7 


1.5 


0.11 


0.09 


-2.0 


450 


62 


0.11 


0.38 


0.2 


0.09 


0.06 


0.0 


0.0 


-2.35 


90 


48 


1.2 


16 


5.5 


11 


7.9 


0.47 


0.37 


-2.35 


190 


41 


0.55 


4.0 


1.7 


1.6 


1.3 


0.07 


0.07 


-2.35 


450 


35 


0.04 


0.2 


0.09 


0.07 


0.05 


0.0 


0.0 



Note. — R 6 formation and merger rates for BH-BH , BH-NS , NS-NS and BH-PSR systems, a is the index of the 
IMF power-law, a provides the dispersion of the Maxwellian distribution of kick speeds, and the kick cutoff gives 
the maximum permitted kick speed (acting on a 1.4 Mq neutron star). Black hole kicks are permitted below a 
ZAMS mass of 40 M , and are scaled to the compact object mass. The BH-PSR column refers to a subset of the 
BH-NS systems, where the pulsar is recycled by mass transfer from the black hole progenitor star. The mass transfer 
should in many cases spin-up the pulsar to millisecond-order periods. 



Table 2 



a 


a (km s 1 ) 


BH-BH 


T < T h 


BH-NS 


T < T h 


NS-NS 


T < T h 


BH-PSR 2 


T < T h 


-2.0 


90 


250 


0.81 


31 


9.0 


10 


7.8 


0.63 


0.53 


-2.0 


190 


240 


0.62 


12 


5.5 


1.8 


1.5 


0.08 


0.06 


-2.0 


450 


240 


0.77 


2.0 


1.1 


0.04 


0.03 


0.0 


0.0 


-2.35 


90 


170 


0.55 


21 


6.5 


11 


8.1 


0.55 


0.47 


-2.35 


190 


170 


0.71 


7.8 


3.4 


1.7 


1.4 


0.07 


0.07 


-2.35 


450 


170 


0.32 


1.4 


0.83 


0.04 


0.04 


0.0 


0.0 



Note. — Rq formation and merger rates. Identical to table 1, save that black hole kicks are not permitted under any 
circumstances. This clearly elevates the number of BH-BH and BH-NS systems seen. Interestingly, the application of 
natal kicks increases the fractional merger rate in BH-BH systems, though the large number of BH-BH systems which 
remain bound when kicks are not present means that the number of merging systems is still comparable, especially 
for the most energetic kick parameters. 
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12 13 14 15 16 12 14 16 



semi-major axis (log cm) semi-major axis (log cm) 

Fig. 1. — A plot of the semi-major axis versus ec- 
centricity for all BH-BH (left) and BH-NS (right) 
systems generated in a single run of 10 6 binaries. 
Triangles designate systems which will coalesce in 
less than 10 Gyr, whilst circles are systems sta- 
ble against gravitational coalescence. For the BH- 
BH systems, a clear distinction between merging 
and stable binaries is evident. The bulk of those 
BH-NS binaries which remain bound after the sec- 
ond supernova also merge within a Hubble time. 
The IMF index used is a = 2.35, the Salpeter 
IMF. The dispersion in the Maxwellian kick ap- 
plied at each supernova is a = 90 km s _1 , with 
the kick maximum K rnax = 500 km s _1 . All kick 
speeds apply to a 1.4 M Q neutron star; the actual 
Aw is found by scaling the kick linearly with the 
compact object mass. 
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BH mass (M.) BH mass <M a ) 



Fig. 2. — Two histograms displaying mass distri- 
butions for BH-BH pairs that remain bound after 
the second supernova. The parameters of the run 
shown are an IMF parameter a = 2.35, kick pa- 
rameters o — 90 km s _1 and K max = 500 km s _1 . 
The left histogram shows that black hole masses 
clump in two locations; one around 5 M Q , and 
the other at 35 M©. We associate the lower- mass 
BHs with 20-40 M Q ZAMS progenitors, which 
lose much of their mass through a supernova (and 
hence experience a natal kick). The high- mass 
set came from more massive progenitors, and col- 
lapsed to a black hole without a supernova or the 
resulting kick. 
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chirp mass (M @ ) 

Fig. 3. — Shows the distribution of chirp masses 
for all BH-BH systems that will eventually merge. 

5 /6 

The gravitational wave amplitude scales as M c ^ irp , 
so this parallels the distribution of signal strengths 
one could expect from a BH-BH population. 




Fig. 4. — A plot of merger time versus the final 
velocity of the BH-BH binaries (left), and BH- 
NS binaries (right). The natal kick parameters 
were a — 90 km s^ 1 , K max = 500 km s _1 (upper 
panels), and a = 190 km s _1 , K max = 500 km 
s _1 (lower panels). Of the systems that remained 
bound after the second supernova, there is a clear 
correlation between increasing system speed and 
a shortening of the gravitational radiation merger 
timescale, from kick-induced high eccentricity or- 
bits, with enhanced GW radiation at periastron 
passage. 
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